Hydrogen atom in phase space. The Kirkwood-Rihaczek representation 
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We present a phase-space representation of the hydrogen atom using the Kirkwood-Rikaczek 
distribution function. This distribution allows us to obtain analytical results, which is quite unique 
because an exact analytical form of the Wigner functions corresponding to the atom states is not 
known. We show how the Kirkwood-Rihaczek distribution reflects properties of the hydrogen atom 
wave functions in position and momentum representations. 
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I. INTRODUCTION 



' The Wigner function [Q, which was introduced in 1932, is up to now the most commonly used phase space quasi- 
distribution representation of physical systems in momentum and position representations. This function fully char- 
^ , acterizes the quantum state and gives basic physical intuition about the investigated state. Moreover, there are 
' some simple physical properties that make this function unique Nevertheless, over the years many others quasi- 
■ distributions have been introduced and studied. Especially well-known are functions from the Glauber and Cahill 
i: s-parameterized class of distributions |^ containing the Wigner function, the Glauber-Sudarshan P-representation 
y—^ |§, H , and the Husimi or the Q-representation |H . 

For many systems the Wigner function is well known and has been carefully studied, but there are some elementary 
quantum mechanical states, for which this phase space representation is not known. An example of such system is the 
hydrogen atom. No exact analytical formula of the Wigner function is known even for the Is state. In the literature 
one can find a limited number of papers devoted to this subject |0, ||, ||, 0, and all the published results have been 
^ ' achieved using methods of approximation. 

^1 It is the purpose of this short communication to present an analytical phase space representation of the hydrogen 
qh, atom using a lesser- known Kirkwood-Rihaczek (K-R) distribution function ||ri| . Its form allows us to obtain analytical 
^jj^ ' results, and as we have shown in [|l^ this distribution, like the Wigner function, fully characterizes a quantum state. 

^: 

. II. DIFFERENT QUASI-DISTRIBUTION AND OPERATOR ORDERING 

In general, the problem of defining a joint probability function in phase space lies in the operator ordering. Because 
different orderings of position q and momentum p operators are not equivalent, the association of the quantum 
operator with classical functions has no unique solution. The association of classical phase space functions A{q,p) 
with quantum operators is shown by: 



{-4(9,p)}ordcring = // dqdp A{q,p) {5{q - q)5{p - p)}„d. 



cnng- 



The Fourier decomposition of the Dirac delta functions makes it possible to rewrite the above equation in the form 
{S{q-q)S{p~p)} ordering^ // dqdpe ^ \e~e ~\ . (1) 

[ZTrtl)'^ J J (. J ordering 

For every operator ordering there exists a corresponding probability quasi-distribution in phase space. 



2 



The Wigner distribution function is obtained by taking a quantum average of the above formula with the respect 
to the Wigner- Weyl ordering: 

W{q,p) = (* I {5{q-mv-p)}ww I *) = JJdq'dp'e-^''^{^ \ e"^^^ I (2) 

in which q and p operators are put in the same exponent. 

The K-R distribution function corresponds to an ordering caUed the anti-standard ordering, that is obtained by 
putting aU p operators on the left of all q operators: 

K{q,p) = I 6{p~mQ^q) I *) = J J dq'dp'e'^'^i^ I e^e-^ | (3) 

The complex conjugation of the K -R function corresponds to the standard ordering (all operators q are on the left 
followed by p operators) , and the real part of the K-R distribution function is obtained by a symmetric superposition 
of the anti-standard and the standard ordering: +p™q"). 

Obviously, the Wigner- Weyl ordering and the anti-standard ordering are not equivalent, so they lead to different 
quasi-distributions. But both of these orderings are intuitive and simple ones, and appropriate phase space distribution 
functions have the same marginal properties. 



III. THE WIGNER AND THE KIRKWOOD-RIHACZEK DISTRIBUTION FUNCTIONS 

The definition of the famous Wigner quasi-distribution function was firstly presented as: 

VF*(q,p) = ^ y**(g + ^/2)e^ *(q-e/2)de- (4) 

Although for most one dimensional systems formula allows to evaluate the Wigner function relatively easily, in 
three dimensions this problem is more complicated. Integrals become quite cumbersome and sometimes impossible to 
handle analytically. As we have mentioned before, an analytical formula for the Wigner function is not known even 
for Is state of hydrogen atom. 

That is why, in order to investigate the hydrogen atom in phase space we shall use the K-R distribution. This 
fmiction was introduced by Kirkwood just one year after Wigner introduced his function. Then, in 1968, the 
same function was rediscovered by Rihaczek [^3| in the context of signal time-frequency distributions and is known 
by his name in signal transmitting theory. The definition of the K-R function in terms of anti-standard ordering has 
a simple form: 



K{q,p) = ^/ de**(e)e'^*(^) 



= ^^iq)e-'^^*{p). (5) 

Elsewhere [ p^ we have presented an extensive analysis of the K-R distribution and its comparison with the Wigner 
distribution function. Here we will only mention the main properties of the K-R function. Similar to the Wigner 
function, the K-R distribution has the correct marginal properties: 



J Kiq,p)dp = \^{q)\^ , 



K{q,p)dq = ^\^ip)W (6) 

which is a fundamental requirement for a joint distribution in phase space. Unlike the Wigner distribution, the K-R 
function is a complex function, and its operator form is not hermitian. Nevertheless, for an arbitrary quantum state 
the knowledge of the K-R distribution at every point of the phase space allows a full reconstruction of this state (see 

M). 
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IV. HYDROGEN ATOM IN PHASE SPACE 



The definition of the K-R distribution (j^) allows one to evaluate analytical formulas of this distribution for many 
complicated quantum mechanical states, for which the analytical expression of the Wigner function is not accessible. 

In the following we will present analytical results obtained for states of the hydrogen atom. The wave functions of 
a non-relativistic hydrogen atom are well known, both in position and momentum representation. The Schrodinger 
equation with Coulomb potential can be separated in the spherical polar coordinates (see, e.g. 10). The solutions of 
angular part of the equation are given by the spherical harmonics Yim{9, (p), and of the radial part can be expressed 
in terms of Laguerre functions Lj^^ (in the position representation): 



Rni{r) = - 



(n-;- 1)! 



{n + iy.3{2n) 



2Z\' 
n J 



t21 + 1 



or Gegenbauer functions (in the momentum representation): 



Fnlip) 



2 {n 

TT 



;-i)! 



[n + iy. 



l22(i+l);j 



W+1 

(nV + 1)'+^ "-'"-'--I 



^2p2 _ 2 



1 



(7) 



(8) 



As usual, n, Z, m denote principal, orbital, and magnetic quantum numbers, respectively. 

Substituting the expressions (^ ||) into the definition of the K-R distribution Eq. (H), and inserting the angular 
relation, we obtain the general formulas for K-R functions of the hydrogen atom states. In the next Section we shall 
present graphical results obtained for selected hydrogen quantum numbers. It is worth noting that the absolute square 
of the K-R function: 

1 



{2n) 



is proportional to the product of probability densities in position and momentum representations (h = 1). Thus, one 
can treat it as the cornerstone of the phase space analysis of physical systems. 



V. EXAMPLES 

Figures 1 to 5 present phase space K-R representation of Is, 2s, 2p, 10m, and 10/ states of the hydrogen atom. 
Let us note that for a 3-dimensional systems the K-R distribution is a function of 6 variables {r,9,ip,p,9' ,(p'), so a 
graphical phase space representation can only present selected cross-sections. We have decided to use the following 
method to represent the hydrogen atom in phase space: all the Figures labelled by (a), (b) (later called Fig. X.a, X.b) 
are auxiliary ones - they present hydrogen atom wave functions in position and momentum representation and the 
arrows depicted on those figures show the directions chosen for the cross-sections. The corresponding cross-sections 
of the K-R phase space functions are plotted in Figures labelled by (c) and (d) (Fig. X.c, X.d). 

All the Figures are organized as follows: Fig X.a shows cross-section tp — const of the absolute square of hydrogen 
atom wave function in the position representation, i.e. | Rni{i')Yim{9 , p) P which is the probability density of finding 
an electron at point (r, 9, ip) (and which does not depend on the value of ip). Fig. X.b shows cross-section p' = const of 
the absolute square of the same wave function in the momentum representation, i.e. | Fni{p)Yi„i{9' , p') p which is the 
probability density of finding an electron with momentum (p, 9' , p') (and here again, a cross-section does not depend 
on the value of p'). The arrows depicted on the Fig. X.a, X.b show the directions chosen for the cross-sections of the 
corresponding K-R phase space functions that are presented in Fig. X.c, where the real part of the K-R distribution 
and appropriate contour plot arc shown, and in Fig X.d, where the absolute value of this K~R distribution and its 
contour plot are plotted. 

Let us first study the K-R function of the Is state. Well-known solutions of the radial part of the Schrodinger 
equation are given by, in atomic unit: 

i?io(r) = 2e~^ (in the position representation), and 

Fio{p) = V f (14* '^Y^ i^^ momentum representation). 
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Fig. l.a shows spatial probability density | i?io(r)/(47r) p of finding an electron at certain point, and Fig. l.b shows 
momentum probability density | Fiq(p) / (Air) p of finding an electron with definite momentum. Note, that those are 
not radial distributions and they have nonzero values at r = and p = as the wave functions were not multiplied 
by or factors. The K-R distribution of Is state is given by: 

Kio{r,e,^,p,e',^') - (27r3)-i^^^-^exp(-iprcose), (9) 

where Q denotes the angle between r and p: 

cose = cos{e-9') + {cos{(p-ip') - l) sin 6* sin 6*'. (10) 
Fig. l.c shows the real part of Eq. (||), K^^ is, multiplied by (27r)'^: 



Kne 1. = (27r)^Re 



71 TT 



(11) 



The cross-section is chosen in the directions depicted by arrows in Fig. l.a, l.b. Next we have displayed its contour 
plot. The absolute value of Eq. (||) multiplied by {2tt)^ is: 



KAbs Is = (27r)3Abs 



TT TT 

Kw[ r, ■2>0,p, -,0 



(12) 



and its contour plot is presented in Fig l.d. In this case, K^ic and KAhs look quite similar mainly because the 
wave function both in position and momentum representation decreases rapidly with the increase of r or p, and the 
oscillating cos(rp)-like structure characteristic for the real part of K-R distribution is not clearly seen on the scale 
used to make these figures. These oscillations are merely marked in the contour plot of K^c by dashed lines which 
correspond to i^Rc = 0. 

Next we shall examine the K-R distribution of 2s and 2p states. Solutions of the radial part of the Schrodinger 
equation in momentum representation are given by: 

32 4p2 _ 1 

^20(P) = 



y^(l+4p2)2^ 

128 p 



^"^^^ ~ V3^(l + 4p2)2- 
Thus, the corresponding K-R distribution are of the form: 

4p2 



K2o{r,e,^,p,e',ip') =/(25^(2-r)e-i^^^^-5^exp(-jprcose), (13) 

K2i{r,e,ip,p,0',ip') ^ y-^^^re^i- — I— — exp(~ ipr cos Q) cos cos 0' , (14) 

3 (1-1- 4p^)^ 

where <d is defined by Eq. ([l0|). The real part of Eq. ( |l3|) and appropriate contour plot are shown in Fig. 2.c and 
the real part of Eq. (|lj) in Fig. 3.c. The absolute value of Eqs. (|l^,|lj) and their contour plots are presented in Fig. 
2.d and Fig. 3.d, respectively. The plot presented in Fig. 3.d has one maximum at (r,p) = (2, \/3/6) the location of 
which does not depend on variables 0, ip, 0' , ip' chosen for the cross section. This result is in perfect agreement with 
what we expect for the state which has quantum numbers n, I satisfying the condition n — I = 1. In Fig. 3.c, where 
the real part of the K-R distribution is presented, we see the additional minimum. It is due to cos(rp)-like oscillations 
of the real part of the K-R distribution. Actually, there is an infinite number of alternate maxima and minima but 
their amplitude decreases rapidly and only two are seen on this scale. Dashed lines depicted on this (and every other 
contour-plot) denote Kjje = 0. According to the results of Fig. 2.d the absolute value of the K-R distribution for 
the 2s state also has quite a simple form: we notice one global and three local maxima on this cross section. Fig. 
2.C, where the real part of the same function is presented, is much more complicated when we look at it close up (see 
contour plot): Extrema are localized on lines where cos(rp) achieves its extremal value, but on this regular structure 
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is another one - the wave function dependence on r and p that is responsible for the appearance of additional zeros 
and sign-changes of the plotted figure. 

Finally we shall present an instructive example of a Rydberg state. Figure 13 shows the K-R distribution for a state 
with n = 10 and I = 9. The real part of the K-R distribution is shown in Fig. IJ.c, and the absolute value is shown in 
Fig. |4|d. Just as one would expect there is one peak in the absolute value plot and similar, but modified by cos(rp), 
structure in the real part of the plot. Comparing the figures with that shown in Fig. ^.c-d we notice similarities 
in their structure. Obviously, peaks are located at different values of {r,p) but there is the property that for states 
with n — I ~ I we observe one maximum in the absolute value plot. Its location changes with the increase of n (it 
moves towards larger r and smaller p) and the maximum value gets much smaller, but it is apparent this class of 
states has much in common. Figure |^ shows the K-R distribution for state n = 10, ^ = 8. Analysis similar to that 
made above holds true when we compare Fig. || and ^ Without counting that all the scales have changed we find 4 
extrema ((n — l)'^) in the absolute value plots and similar oscillating structures in the real part plots. It is a general 
and well known result that n — I is equal to the number of extrema of hydrogen atom wave functions both in position 
and momentum representation. We only stress that the K-R distribution renders very well this property: its absolute 
value has {n — l)'^ extrema. 



VI. SUMMARY REMARKS 

We have presented a phase space representation of the hydrogen atom using the K-R distribution. An advantage of 
such a representation is that the system can be treated analytically. It is an attractive feature of the K-R distribution 
that as long as one knows the wave function of the system there is no need to perform any integrals to obtain this 
distribution value in every point of phase space. 

We have described in detail the Is, 2s, 2p, 10m, and lOl states of the hydrogen atom. Plots corresponding to the 
K-R functions of these states were presented and compared with the plots of probability densities in position and 
momentum representations. 
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FIG. 1: The K-R phase space representation of Is state (n = 1, I = 0): 

(a) the cross section ip = const of spacial probability density | _Ri()('r)/(47r) f of finding electron at point {r,6,ifi); 

(b) the cross section (p' = const of momentum probabihty density | Fio(p)/(47r) |^ of finding electron with momentum (p, 0' , (p'); 

(c) the real part of the K-R distribution and its contour plot (the cross section is made in directions shown by arrows depicted 
on Fig. 1(a) and (b)); 

(d) the absolute value of the K-R distribution and its contour plot (the same cross section as above) . 



FIG. 2: The K-R phase space representation of 2s state {n = 2, I = 0): 

(a) the cross section = const of spacial probability density | J?2o(f')/(47r) |^ of finding electron at point {r,0,(p); 

(b) the cross section p' — const of momentum probability density | F2o{p) / (47r) p of finding electron with momentum {p, 9' ,ip'); 

(c) the real part of the K-R distribution and its contour plot (the cross section is made in directions shown by arrows depicted 
on Fig. 1(a) and (b)); 

(d) the absolute value of the K-R distribution and its contour plot (the same cross section as above) . 



FIG. 4: The K-R phase space representation of 10m state (n = 10, Z = 9): 

(a) the cross section ip = const of spacial probability density of finding electron at point {r,9,ip); 

(b) the cross section (p' = const of momentum probability density of finding electron with momentum (p, 6' , ip'); 

(c) the real part of the K-R distribution and its contour plot (the cross section is made in directions shown by arrows depicted 
on Fig. 1(a) and (b)); 

(d) the absolute value of the K-R distribution and its contour plot (the same cross section as above) . 



FIG. 5: The K-R phase space representation of 10/ state (n = 10, Z = 8): 

(a) the cross section tp = const of spacial probability density of finding electron at point (r, 6, ip); 

(b) the cross section ip' = const of momentum probability density of finding electron with momentum {p, 9' , ip'); 

(c) the real part of the K-R distribution and its contour plot (the cross section is made in directions shown by arrows depicted 
on Fig. 1(a) and (b)); 

(d) the absolute value of the K-R distribution and its contour plot (the same cross section as above) . 



FIG. 3: The K-R phase space representation of 2p state ( n = 2, Z = 1 (a) the cross section ip = const of spacial probability 
density | ii2i(r)/(47r) ^ of finding electron at point {r,9,(p); 

(b) the cross section tp' = const of momentum probability density | F21 (p) /(47r) |^ of finding electron with momentum (p, 9' , tp'); 

(c) the real part of the K-R distribution and its contour plot (the cross section is made in directions shown by arrows depicted 
on Fig. 1(a) and (b)); 

(d) the absolute value of the K-R distribution and its contour plot (the same cross section as above) . 
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